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Abstract 

Built on a recent work on the quantum chromodynamic (QCD) factorization for semi-inclusive 
deep-inelastic scattering (DIS), we present a set of factorization formulas for the spin-dependent 
DIS and Drell-Yan cross sections at low transverse momentum. The result can be used to extract 
transverse-momentum dependent parton distribution and fragmentation functions from relevant 
experimental data. 
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1. Polarized hard scattering has becoming an important tool to learn about the internal 
structure of hadrons. In recent years, inclusive and semi-inclusive polarized deep-inelastic 
scattering (DIS) has helped to unravel the quark helicity distributions in the nucleon P|. 
At present the polarized Relativistic Heavy Ion Collider at Brookhaven National Labora- 
tory is producing data from polarized proton-proton collisions, which will provide a direct 
measurement of the polarized gluon distribution 

An important class of polarized experiments involves measurement of transverse momen- 
tum of the order Aqcd- For example, in semi-inclusive deep-inelastic scattering (SIDIS) at 
moderate energy (for example, at HERMES kinematics), the tagged final-state hadron has a 
transverse-momentum peaked at a few hundred MeV. In Drell-Yan process also at moderate 
center-of-mass energy, Drell-Yan pairs typically have a transverse momentum of the same 
order of magnitude. Theoretical study of these processes began with the classical work of 
Collins and Soper in which a nearly back-to-back hadron pair is produced in e+e~ collisions 
A factorization theorem for the process was established, which involves a new class of 
non-perturbative hadronic observables depending on the transverse-momentum of hadrons 
and/or partons: the transverse-momentum dependent (TMD) fragmentation functions and 
parton distributions. In a previous publication 0|, we have extended the QCD factorization 
theorem to the case of the semi-inclusive DIS. The correction to the factorization is on the 
order of P\IQ^ and Mf^jQ^ ^ where Pi is the transverse momentum of the produced hadron 
and Mu is a hadron mass scale. In this paper, we extend factorization further to a class 
of spin-dependent DIS and Drell-Yan processes. Since the details are similar, here we focus 
mostly on the final result, omitting technicalities which can be found in the above refer- 
ences. The result can be used to extract a new class of TMD parton functions from relevant 
experimental data, just like the standard QCD factorization theorems allow extraction of 
the Feynman parton distributions from hard scattering data. 

Before proceeding future, let us remark that there have been many studies in the literature 
on the same processes at large, but not too large transverse-momentum (Q ^ -Pi ^ Aqcd), 
where Q is hard-collision scale (for example, the virtual-photon mass). In this case, the cross 
section can in principle be calculated using the conventional perturbative QCD method with 
the integrated Feynman parton distributions. However, the hard part contains the large 
double logarithms of the type a^ln^Pj/Q^, which must be re-summed to make reliable 
predictions The formalism developed by Collins and Soper in the case of e+e~ annihi- 

lation PI and followed in our previous paper, is ideal for making this type of re-summation. 
In fact, an application to the unpolarized Drell-Yan process was first made by Collins, Soper 
and Sterman (CSS) (2]. Various applications of CSS formalism to Drell-Yan, heavy bo- 
sonproductionS;^ 

Bi mini HQ. 

2. Let us first define a new class of non-perturbative hadronic matrix elements, the 
spin and transverse-momentum dependent parton distributions and fragmentation functions, 
which one hopes to learn from high-energy scattering. As illustrated in Fig. 1, we consider 
a hadron with momentum moving in the z-direction, and is polarized with a spin vector 
5''^ (dimensionless and S" ■ P = 0). In the limit oo, the P^ is proportional to the light- 
cone vector j)^ = A(1,0,0, 1), where A is a mass dimension-1 parameter. The conjugation 
light-cone vector is n = (1/2A)(1,0, — 1), such that = n"^ = and p ■ n = 1. We use 
the light-cone coordinates fc^ = (k^ ± k^)/\^, and write any four-vector in the form 
of {k~,k) = {k~ , k~^ , k±), where k± represents two perpendicular components {k^^k^). Let 
(xP^, k^) represent the momentum of a parton (quark or gluon) in the hadron as shown in 



jroductions, and semi-inclusive DIS scattering have been developed in the literature 
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FIG. 1: Transverse momentum dependent quark distributions in the nucleon. 

Fig. 1. In a non-singular gauge (e.g. Feynman gauge), the TMD parton distributions can 
be defined through the following density matrix [ij, , 



M^{x,k±,n,xC,p) = 



27r 



(27r 



(1) 



X 



PS 



PS 



S^b^,fi^,p) 



where ipg is the quark field, and the Dirac- and color indices of the quark fields are implicit. 
The +(— ) superscript is appropriate for DIS (Drell-Yan) process jisl . EgI ]. v'^ is a time- 
like dimensionless (f^ > 0) four-vector with zero transverse components {v~,v~^,0) and 
v~ ^ v^. Thus the is a quasi light-cone vector, approaching n'^. The variable ("^ denotes 
the combination (2P ■ v)"^ /v"^ = C^. is a gauge link along v^, 



(2) 



Here the non-light-like gauge link is introduced to regulate the light-cone singularities. We 
avoid the use of singular gauges (e.g. the light-cone gauge) because in those gauges the 
gauge potential may not vanish at infinity and gauge links at infinity might be necessary to 
define gauge-invariant parton distributions [l6| . 

In the above definition, we have derived a soft factor defined as 



S^{b^,fi\p) = ^(0|4,(6x,-oo)4^.(±oo;6x)/:,,fc(±oo;0)/:,fc,(0;-oo)|0) , (3) 

where i,j, k, I = 1,2,3 are color indices and new quasi light-cone vector v'^ = {v^, v^, 0) has 
been introduced with v~ <^ v^. The p parameter is defined as p = -^Z v~v^ /v^v~ 3> 1. The 
above soft factor will also be present in the factorization theorems below. 

The leading order expansion of the density matrix M. contains eight quark distributions 



17, 1 



M 



q{x, ki_) ^ + j^^Qi^, k±)a^''ki,p^ + Agi(x, k±)X-f5 i> 



1 - - 1 

^—l^qrix, k±)-f5 ^(A;_L ■ S±) + —SgL^x, k^^Xia^^-f^p^k^i + 5qT{x, k±)icrf,^-f5p>'S'i 



+ ^^SqT'{x, k^)ta^,^,p^ (k^ ■ S^k'i - IklSl] + ^qrix, k^)e^-"'''^^p,k^Sp 



M2 



(4) 



where M is the nucleon mass. We have omitted ± labels and the arguments C, and /i for the 
distributions at the right side of the equation. The convention for 75 and e-tensor follows that 
of [19|J. In principle, there are Lorentz structures depending on and they can taken 
to be and p^, respectively, when no light-cone singularity is present. The polarization 
vector has been decomposed into a longitudinal component Si and a transverse one S"^, 



where A is the helicity. The notations for the distributions follow Ref. |12(j], which are 
different than those in 

Out of the eight TMD distributions, three of them are associated with the k±-even struc- 
ture under the exchange kj_ —k±: q{x, kj_), Aqii^x, k±), and Sqrix, k±), which correspond 
to the unpolarized, longitudinal polarized, and transversity distributions, respectively j2ll |. 
These distributions will survive after integrating over transverse momentum. The other 
five distributions are associated with the k±-odd structures, and hence vanish when k± 
are integrated. Two of them, qT(yX,k±) and 6q{x,k±), are odd under naive time-reversal 
transformation ^23\ and require hadronic final-state interactions to be non-zero @]. 

If one is interested in production of hadrons, TMD fragmentation functions must be intro- 
duced. Following the above, we define the following density matrix for a quark fragmentation 
into a (pseudo) scalar hadron, 

M,{.,P,.,,,(/z.p) = !f / |-^e-"=*^-«-'-J (6) 

X a 

X (p,x| (v^„,(r, 6)4(r, b, -00) |o)/5(6^, fi, p) , 

where v is mainly along the light-cone direction conjugating to Ph, k~^ = / z and = 
—Ph±/z, and a = 1, 2, 3 is a color index. The variable C, is defined as = 4(P/j ■ /v'^. At 
leading order, we can have two fragmentation functions in the expansion. 



g(x,p_L) 7^- + j^Sq{x,p±)a^"'pf,±n^ 



(7) 



The second one is (naively) time-reversal odd. 

The fc^-even parton distributions satisfy the same CoUins-Soper evolution equation re- 
spect to C, e.g., jsf 

d 

C^/(a;, b, /i, x() = {K{h, n) + G{x(, fi))f{x, b, p, xC) , (8) 

where f = q, AqL, and 6qT, and the sum K + G is independent of ultraviolet scale /i. This 
is because the ^-dependence arises from light-cone divergences that are independent of the 
spin structure. We have verified this explicitly at one-loop order. For k±-odd distributions, 
similar equations have not yet been studied in the literature. 
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3. In j^, a factorization formula has been shown for the cross section of unpolarized semi- 
inclusive DIS, which involves the unpolarized TMD parton distributions and fragmentation 
functions. Let us first extend this result to polarized semi-inclusive DIS process. Our result is 
valid up to power corrections in P]_/Q'^, and can be used to extract the new non-perturbative 
hadronic matrix elements from relevant experimental data. 

The differential cross section for semi-inclusive DIS reads 

-^^"'-l/£,.l^-^(P,g,P.), (9) 



dxBdydzhCpPhi. Q 
where the leptonic tensor is 

= 2{Fr + tr - ^^"gV2 - 2i\te^""'^i'jp) , (lo) 

and is the helicity of the initial lepton, £ and ^' the initial and final momenta of the lepton, 
= i'^—i'^ the momentum of the virtual photon, Ph the momentum of the observed hadron. 
As usual, xb is the Bjorken variable, y is the fraction of the lepton energy loss, and = —q^. 
The hadron tensor has the following expression in QCD, 

W^^{P,q,Pn) = ^E / (0ie^'-«(P5|J,(O|XP,)(XP,|J.(O)|P5) , (11) 

where is the electromagnetic current of the quarks, X represents all other final-state 
hadrons other than the observed particle h. The variable Zh can be defined as P ■ Ph/ P ■ q 
or P/T/?"- 

In a coordinate system in which the virtual-photon and hadron momenta are collinear 
and along the ^-direction, the above hadronic tensor W^^ has the following leading-twist 
structures, 

+t\^e'i^i\FLL + S^PI^Flt) 
j-fnt^'^ _ _ p*^ ^ f(^) 

+^Ph±a^'±^ Ph]_Fu L 
nf^'^^cxPc p 171(1) 

+Ph±a^'±^ Ph±{Ph± ■ S±)FI)^ , (12) 

where g'^ = g^^ — p^n^ — p^n" and e^^ = e^^^^^Panp, and P/i^ is the unit vector along P^^. 
Fi^j^ are structure functions depending on x^, Zh, Ph±, and Q^, where Ji represents the 
polarization of the incident lepton and I2 that of the target hadron {U : unpolarized, L: 
longitudinally-polarized, T: transversely-polarized). For example, Fm; denotes the unpo- 
larized structure function. Fur means the target is transversely polarized while the beam 
lepton is unpolarized. 
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Substituting the hadronic tensor into the differential cross section, we get 

, /f = ^^[(l-y + yV2)a:B(F«+sin(</.,-05)|5jF«) 

+Xey{l - y/2)xB{XFLL + cos{<Ph - <Ps)\S±\Flt) 
+ (1 - y)xB (- cos(20OFi'j} + A sm{24>h)FuL 



(13) 



where s = {P+£)^, and (ps and (ph are azimuthal angles for the transverse polarization vector 
of the initial hadron and the transverse momentum of the final state hadron, respectively. By 
studying the angular dependence of the cross section, one can isolate the different structure 
functions. 

Following we can factorize the structure functions into TMD parton distributions 
and fragmentation functions, and soft and hard parts. For example, the double polarized 
structure function Fll has the following factorization form. 



q=u,d,s,... 



X 



Hll {(flp\ p) 5\zk^ +p^ + £l - PhA_) . (14) 



where we have chosen a coordinate system in which = {Q'^/x^)p and ("^ = {Q'^zf^p. 
Similarly, we can have for the other structure functions, 

Flt = J '^^^AqT{xB,k^)q{zf,,p^)S^{i^)HLT{Q') , 
4't = J ^^^qT{xB,k^)qizH,P±)S^{i^)Hi;l{Q') , 

f!j^ = I ^^j^5qT{xB,k^)5q{zH,P±)S^{i^)H^Jl{Q') , 
(3) f 4{k± ■ Ph±fp± ■ Ph± - \k±\^P± ■ Ph± - 2k^ ■ p±k^ ■ Ph± 



/ Mm, 



h 



x6q'T{xB,k^)6q{zh,P±) S''{i^)Hgl{Q') , (15) 

where the simple integral symbol represents a complicated integral used above, i.e., J = 
J d'^k±d'^p±d'^i±5'^{zk± + p± + i± — Ph±)- And for simplicity, we also omit the explicit sum 
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over all flavors weighted with their charge square, and the explicit /i, p and C, dependence for 
the parton distributions and fragmentation functions, as well as the soft and hard factors. 

From the above factorization formulas in the transverse momentum space, we can also 
obtain the factorization form in the impact parameter space. For example. 



■1 



= Jj^^di [{dlqTixB, ZHb)) q {z„ h) S^{h)H)j), {Q^)\ , (16) 

where Fl)^{b) = J d'^Ph±e^^''^'''\Ph±\Fl)'l{Ph±), and other quantities depending on 6 are ob- 
tained from the Fourier transformation of the corresponding momentum-dependent ones. 
The convolution in transverse momentum space becomes a simple product in the impact 
parameter space, while the k± moment integral generates a derivative on b. 

An explicit one-loop calculation for Fuu and Fll can be done. The hard parts at one-loop 
order have the following form. 



(1 + lnp2) In^ - lnp2 + 1 In^ p^ + yr^ - 4 

^ ' jJL^ 4 



(17) 



To get the next-to-leading order corrections for the hard parts related to fc^-odd distri- 
butions, one has to consider two-loop calculations for scattering off an elementary quark 
target. 

In the above discussions, we consider the factorization of the leading contributions in 
the expansion of Ph^/Q"^- However, the power corrections in the semi-inclusive DIS 
might have sizable effects at moderate range. They are not included in our factorization 
formalism. On the other hand, it will be interested to extend our factorization formalism to 
include the sub-leading power contributions in SIDIS, just like what have been done for the 
sub- leading power correction to the inclusive DIS and DY cross sections ^]. This however 
is beyond the scope of the present paper. 

One can apply the above formalism to study the polarized cross sections and asymme- 
tries in the SIDIS, and compare theoretical predictions with experimental measurements if 
the non-perturbative parton functions are known from solving non-perturbative QCD. In 
practice, we can treat them as unknown inputs, and flt to the experimental data. The par- 
ton functions determined phenomenologically can be used to make predictions for similar 
processes. Note that in the TMD quark distributions and fragmentation functions there 
is a double logarithmic dependence on the hadron energy (in terms of In^ &^C^) which is 
controlled by the Collins-Soper evolution equations as shown in Eq. (jSj). To get the reliable 
prediction for the cross sections for the SIDIS, one has to solve these evolution equations 
to re-sum the double logarithms. As we showed in the above the Collins-Soper evolutions 
for the fc_L-even TMD quark distributions are the same. It will be interested to study the 
Collins-Soper evolution for all the leading-twist quark distributions including the /c_L-odd 
ones. We leave this study in a separate publication. 

Let us flnally remark on the two special cases related to the above results. The first case 
concerns cross sections integrated over the transverse momentum of the hadrons Since 
the above factorization formulas are valid only when P^x -C Q, one cannot use them to 
integrate out all Ph^- To do that, one must have factorization formulas valid at all Pux, 
including all power corrections. Alternatively, one can prove the factorization theorems in 
the P/i_L-integrated form, which can be done, but is beyond the scope of this paper. 

The second case concerns the region of transverse momentum Aqcd Ph± -C Q. In this 
region, the above factorization formula is still valid. However, because Ph± is now hard, all 
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structure functions can be further factorized in terms of ordinary twist-two and twist-three 
parton distributions. Using Colhns and Soper equations, one can sum over large logarithms 
of type a,ln2QVP^\. This case, has been studied for example in Some of the 

relevant leading-order perturbative QCD calculations have been done in j25|. 

4. In the Drell-Yan processes, we have two hadrons in the initial states: A and B. We 
consider the hadron A moving along 2;-direction while B along -z-direction. For the TMD 
parton distributions in hadron B, we introduce another non-light-like vector v = {v~, v'^, 0) 

{v^ <^ v'^) and the variable ( = {2Pb ■ vf' j'xP'^ just like for the fragmentation functions 
in DIS. The TMD quark and anti-quark distributions for hadron B are similar to those of 
hadron A. 

We study the the production of a lepton pair with low transverse momentum in Drell-Yan 
process: 

A(Pi, Si) + i?(P2, ^2) ^ 7*(g) + ^ - ^+ + r + X, (18) 

where P\ and P2; 5*1 and 5*2 are momenta and polarization vectors of A and P, respectively. 
The differential cross sections for t^t~ production reads 

'^'^ "''^.L^.W^'^'^ , (19) 
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where Lp,y is the leptonic tensor. For the Drell-Yan production, we have L^^ = 4(^5*^2 + 
^1-^2 ~ 1'^)) where Ci is the lepton's momentum. VL is the solid angle of the lepton 

in the virtual photon rest frame. is the momentum of the lepton pair; s is the total 
center-of-mass energy square. The hadronic tensor is defined as: 

W^^.(xi,X2,Qx,5i,S2) = j -0^,e-'''-HPiSiP2S2\MOJM\PiSiP2S2) , (20) 

where is the electromagnetic current. Disregarding the contributions from k±-odd parton 
distributions, the leading hadronic tensor can be decomposed into three terms, 

= -g^j_'{Wo - X1X2WLL) - g^T^TT , (21) 

where gl^j^ = Si± ■ Si±g'^ + S'f S'2 + S'(S2- The tensor structure Wq denotes the unpolarized 
tensor, Wll the double longitudinal polarized tensor, and Wtt the double transversely po- 
larized tensor. We can choose a frame that the momenta Pf, P^ and can be decomposed 
into the following forms, 

P^ = Jn'^ + —P' , 
2 s 

Q/^ = Xip'^ + XaV + Q^^ . (22) 
The polarization vectors of hadrons can be decomposed as: 
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where Ai and A2 are the hehcities of hadron A and B respectively, and Si± the relevant 
transverse polarization vectors. 

Using the hadronic tensor, the differential cross section can be expressed as 

-(1 + COS^ e)XiX2WLL{xi,X2, Q^ Q±) 

+ sin2 e cos(0i + <P2)\Si±\\S2±\Wtt{xi, X2, Q\ Qs_) } . (24) 

where Q is the polar angle of the lepton momentum, and 0i and 02 are azimuthal angles of the 
nucleon transverse spin relative to the lepton plane. Wj^i^ contributes to double longitudinal 
spin asymmetry and Wtt to double transverse spin asymmetries. 

Following Ref. 0], we have the following factorization formulas for the above hadronic 
tensors. 



e2 f 

Vro(xi,a;2,Q',Q±) = f d''h^d''k2±d^ 

^ „ ^ 



q=u,d,s^. 



xq{xi, ki±, f/, XiC, p) q (xa, k2±, /i^ X2C, p) S (/l, P^ p) 
xHo{Q\p\p)6\h± + k2±+i±-Q±) , (25) 

WLLixi,X2,Q'',Q±) = Yl f I d''k,^d''k2±d^ii. 



e2 r 

Wtt{xi,X2,Q^Q±) = f I d^h±d^k2±d 



q=u,d,! 

xAqL (xi, kij_,p'^,XiC,p) Aq^ (x2, k2±, p^,X2(, p) S~{i±,p^,p) 
xHll {Q\p\p) S\h± + k2± + L - Q±) , (26) 

— I a Ki±U tv2±u'^' 

q=u,d,s,. 

x5qT (xl,fcl_L,p^a;lC,p) Sqj^ (^2 , A;2± , /i^ X2C, p) 5'"(^L,/i^p) 
xHtt {Q', p\ p) + k2L + /± - Q±) , (27) 

2 

where we have chosen a special coordinate system: = {Q'^ /x\)p and C = {Q'^ /x'^p and 
p = \Jv~v^ /v^v~ . The above results are accurate at leading powers in {Q\/Q'^) for soft 
~ Aqcd- Introducing the impact-parameter space representation, the above factorization 
formulas become simple products, e.g., for the unpolarized tensor, 

wo(xi,x2,6,(5^) = Y (^l'^'/^^^lC,p)^(a;2,&,p^x2c,p) 

q=u,d,s,... 



X 



S~{h,p\p)H,{Q\p\p) . (28) 



The other two hadronic tensors can be obtained similarly. 

In the following, we will show that the above factorization formulas are valid at one-loop 
order, and the relevant hard parts emerge from the calculation. We first consider the soft 
factor in the factorization formulas. Since it has no spin dependence, it contributes equally 
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FIG. 2: One-loop corrections to DY processes. 

to all three structure functions. From its definition Eq. (jSj), the one- loop contribution Drell- 
Yan soft factor is, 



s (^x,/i^p) 



27r2 



In 



4(f ■ V 



1 



+ A2 



nS\i^) In 



(29) 



At this order, it is the same as that for the DIS process 

One-loop Feynman diagrams for the Drell-Yan process initiated from a quark and anti- 
quark pair are shown in Fig. 2, which contain both virtual and real corrections. The virtual 
corrections are the same for all three structure functions, 



Wo,LL,TT = l^ix, - 1)6{X2 - 1)S\Q^){1 + 2{Zf - 1) + 2{Z{. - 1)) 



(30) 



where Zp and Zy come from self-energy and vertex corrections respectively. Zp is the same 
as that in DIS j^. The vertex correction is different, 



a.. 



Air 



Q 



m 



Z{, = l-T^C^ ln-^ + ln'-^ + 21n — ln-^-41n-^- 



47^' 



A2 



(31) 



which can be obtained through analytical continuation from space-like region (g^ < 0) to 
time-like region (g^ > 0). From the above, the one- loop TMD parton distributions, and the 
soft factor in Eq. (jSHl), the one- loop contribution to the hard parts can be extracted. 



(l + ln/)ln^ 



lnp2 + ^In^ p2 ^27r2 -4 
4 



(32) 



where we have chosen a coordinate system in which xiC, = X2C and therefore the dependence 
on the quasi-light-like vectors v and v is simply through a combination p = v~v^ lv^v~ . 
Note that the hard parts at this order are the same for all three structure functions and 
hence are spin-independent. 
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The real contributions are different for the three hadronic tensors: 



_ asCp _ . r 1- xi 2xi(l - xi)m^ 

°~ 67r2 ^ ^^>\q2^ + xiX^ + m^{l-Xi)^ {ql + XiX^ + m^{l - Xi)^^^ 



+ 7-. ^ T^. TT. + ^(1 - Xi) „ , ^ M n- 



^1 — a;i)+ + xiA^ + m2(l — + gl + A^ 

+ (a;i X2) , (33) 

a^Cp r 1 — Xi 2m'^{l — xi){l — xi + xl) 



Gtt^ I g^ + XiA^ + m2(l — (g^ + XiA^ + m2(l — xi~ 

, 2x1 1 , 1 , 



2^2 



l-xi)+gi + XiA2 + m2(l-xi)2 ' ' gi + A2 g^ + A^ 

+ (xi ^ X2) , (34) 
_ ttsCF.,. . r 2xi 1 

DTT^ (1 " Xi)+ gl + 3;iA^ + m^(l — Xij^ 

2xi(l-xi)m2 , 1 , Q2 



(g2 +xiA2 + m2(l-xi)2)2 ^^^^ ^i)g2 ^A2^'' g2 +A2 
+ (xi ^ X2) . (35) 

All these results can be reproduced by one-loop parton distributions for the quark and 
anti-quark and soft factors discussed before. Therefore, the factorization formulas for the 
Drell-Yan processes at low transverse momentum are correct at one-loop order. For an 
all-order proof, one can follow the same procedure outlined in 0]. 

Finally, we comment that when the transverse-momentum of the lepton pair is large 
compared to Aqcd? but much smaller than Q, the above formalism leads to a resummation 
of double logarithms as in the original paper by Collins, Soper, and Sterman m. For spin- 
dependent part, some studies along this direction can be found in Refs. |26l l27|. 
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93ER-40762. J. P.M. was supported by National Natural Science Foundation of P.R. China 
through grand No. 19925520. X. J. is also supported by an Overseas Outstanding Young 
Chinese Scientist grant from NSFC. 
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